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We study the evolution and saturation of the gluon distribution function in the quark-gluon 
plasma as probed by a propagating parton and its effect on the computation of jet quenching or 
transport parameter q. For thermal partons, the saturation scale is found to be proportional to 
the Debye screening mass /if,. For hard probes, evolution at small x — Ql/GET leads to jet energy 
dependence of q. We study this dependence for both a conformal gauge theory in weak and strong 
coupling limit and for (pure gluon) QCD. The energy dependence can be used to extract the shear 
viscosity r] of the medium since can be related to the transport parameter for thermal partons in 
a transport description. We also derive upper bounds on the transport parameter for both energetic 
and thermal partons. The latter leads to a lower bound on shear viscosity-to-entropy density ratio 
which is consistent with the conjectured lower bound r?/s > l/47r. We also discuss the implications 
on the study of jet quenching at the BNL Relativistic Heavy Ion Collider and the CERN Large 
Hadron Collider and the bulk properties of the dense matter. 
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I. INTRODUCTION 

Experimental data from the Relativistic Heavy-ion 
Collider (RHIC) have shown significant suppression of 
both high transverse momentum single inclusive hadron 
spectra and back-to-back dihadron correlation in cen- 
tral high-energy heavy- ion collisions P, [2i The ob- 
served jet quenching phenomena can be attributed to 
parton energy loss and medium modification of the effec- 
tive parton fragmentation functions 0, S 0] due to gluon 
bremsstrahlung induced by multiple parton scattering. 

Within the picture of multiple parton scattering in 
QCD, the energy loss for an energetic parton propagat- 
ing in a dense medium is dominated by induced gluon 
bremsstrahlung. Taking into account of the non-Abelian 
Landau-Pomeranchuck-Midgal (LPM) interference, the 
radiative parton energy loss |7|], 



a,Nr 



(1) 



is found to depend quadratically on the medium length 
L and a jet transport or energy loss parameter. 



I , 2 C^O--R 2 



(2) 



which describes the averaged transverse momentum 
transfer squared per unit distance (or mean-free-path). 
Here R is the color representation of the propagating 
parton in SU{3) and p is the color charge density of 
the medium. According to this picture, jet quenching 
as observed in high-energy heavy-ion collisions is a direct 
measurement of the jet transport parameter in dense 
medium which not only characterizes the color charge 
density but also the interaction strength between the 
propagating parton and the medium. 

Phenomenological studies based on variations of the 
parton energy loss picture d, H, [13, El all indicate the 
formation of an extremely high density matter in the 



initial stage of high-energy heavy-ion collisions at the 
RHIC energy. The averaged transport parameter ex- 
tracted from different phenomenological studies of the 
single inclusive high pT hadron suppression in the most 
central Au + Au collisions at RHIC is g^? ^ 1 — 15 
GeV^/fm (for a propagating quark) at an initial time 
To = 1 fm/c. A recent simultaneous fit of the next- 
to-leading order (NLO) pQCD calculation to both sin- 
gle and back-to-back dihadron suppression ^12] narrows 
the uncertainty to qp — 1-1 ~ 1-4 GeV^/fm, which is 
still about 100 times higher than that in a cold nucleus 
qp ~ 0.013 GeV^/fm as extracted from leading hadron 
suppression in deeply inelastic scattering off large nuclei 

Q. 

In most of the theoretical studies of parton energy 
loss 0, EE EE [13, [3 J except the twist-expansion ap- 
proach [19[, a static potential model for jet interaction 
with the medium was assumed which led to the factor- 
ized dependence of parton energy loss on the transport 
parameter q^j in Eq. ([T]). In this static potential model, 
energy and longitudinal momentum transfer between a 
jet parton and the medium is ignored. Therefore, elastic 
energy loss due to the recoil of the medium parton dur- 
ing the jet-medium interaction is neglected in the cal- 
culation of radiative parton energy loss. Furthermore, 
the static potential model does not include the effect of 
inelastic break-up (or parton radiation) of the medium 
partons which can give rise to jet energy dependence of 
the transport parameter qn. In a dynamical picture, the 
transport parameter can be related to gluon distribution 
density of the medium Q- The jet energy dependence 
of the transport parameter is then directly related to the 
scale and momentum fraction dependence of the gluon 
distribution density. Understanding the jet energy de- 
pendence of the transport parameter not only helps us 
to improve the phenomenological study of experimental 
data on jet quenching but also provides additional infor- 
mation about the structure of the dense quark-gluon mat- 
ter in heavy-ion collisions. Furthermore, as illustrated in 



2 



a recent study [2^ , the low energy limit {E ^ T temper- 
ature of the medium) of the transport parameter in jet 
quenching is directly related to the shear viscosity of the 
quark-gluon matter in a transport description. There- 
fore, experimental and theoretical study of the jet energy 
dependence of the transport parameter will be able to 
provide another piece of important information on bulk 
properties of the dense medium. 

Recently, the transport parameter qn has also been 
calculated for a strongly coupled A/" = 4 supersymmetric 
Yang-Mills (SYM) theory in different limiting scenarios. 
With a definition in terms of an adjoint Wilson loop along 
the light-cone, Liu, Rajagopal and Wiedemann [21|] found 
that QA in the large limit of the t' Hooft coupling A = 
Ncg^ in SYM, 



7r3/2r(3/4) 

r(5/4) 



(3) 



scales with the temperature cubed and is independent of 
the jet (propagating parton) energy. In another limit for 
a slowly moving heavy quark, the transport parameter. 



(4) 



as defined in Eq. ([2]), is found [H, HI], to depend on 
the square-root of the heavy-quark energy, where 7 = 
E/M < (Af/%/Ar)2. It is not clear how these two re- 
sults are related to each other, though both describe the 
transport properties in a SYM theory. 

In this paper, we investigate the jet energy dependence 
of the transport parameter qn within perturbative QCD 
(pQCD). We will first re-exam the relationship between 
the transport parameter and the unintegrated gluon dis- 
tribution function of the color charges in the medium 
and how they are related to parton energy loss in the 
medium. For energetic jet partons, there are large log- 
arithms of both momentum scale and small momentum 
fraction. They allow us to take double logarithmic ap- 
proximation (DLA) and resum gluon radiation of the tar- 
get color charges to all orders. The initial condition to 
such a resummed evolution of the gluon distribution can 
be calculated perturbatively within pQCD at finite tem- 
perature with hard thermal loop (HTL) resummation. 
From such resummed gluon distribution one can further 
take into account gluon saturation and calculate the sat- 
uration scale self-consistently which will determine the 
transport parameter and its jet energy dependence. 



II. PARTON ENERGY LOSS, GLUON 
DISTRIBUTION FUNCTION AND 
TRANSVERSE MOMENTUM BROADENING 

Multiple parton scattering within the high-twist ex- 
pansion framework [IJ] can go beyond the static potential 
model and include energy and longitudinal momentum 
transfer in the calculation of medium modified fragmen- 
tation functions. The total energy loss for a propagating 




FIG. 1: Feymann diagram for induced gluon radiation that 
contibutes to the quark energy loss. 



parton in a deeply inelastic scattering (DIS) off a large 
nucleus as shown in Fig. [1] due to secondary quark-gluon 
scattering in this framework can be expressed as [l9| 



AE 2asCR f d^qr 



dz 



1 + (1-^)2 T^A{xB,XL,qT) 



where 



XL 



2z(l - z)p-k' 



(5) 



(6) 



is the total longitudinal momentum transfer related to 
induced gluon radiation with final transverse momentum 
(.T- A similar longitudinal momentum transfer 



XT 



2qT • It 



2(1 - z)p-k 



(7) 



is always provided by the initial gluon with transverse 
momentum qx- As illustrated in Fig. [11 p = [0,p~,Ot] 
is the initial quark momentum after its interaction with 
the photon, k = [fc"*", 0, Ot] is the momentum per nucleon 
in the medium, qx is the transverse momentum of the 
gluon exchange with the medium, It is the transverse 
momentum and z is the fractional longitudinal momen- 
tum carried by the radiated gluon with four-momentum 
£ = [i'^/2zp~ , zp~ ,£t]- The quark distribution function 
fqixs) represents the production rate of the initial quark 
carrying x b fraction of the nucleon longitudinal momen- 
tum in DIS. Eq. ([5|) is derived for quark energy loss and 
one can extend it for gluon by replacing the correspond- 
ing Casimir factor Cr for gluons. In the coUinear ex- 
pansion of the twist expansion approach, one normally 
makes Taylor expansion of the hard partonic parts in qx 
and only the quadratic terms lead to the twist-four con- 
tribution. One can, however, approximate higher twist 
contributions from the g^-dependence of the hard par- 
tonic part of the multiple scattering by using the average 
value {q^) = fJ-T in the cross section. As an extension 
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of the twist expansion, we will keep the integration over 
the gluon's transverse momentum qx- The unintegrated 
quark-gluon correlation function is defined as, 



Tqgix,XL,qT) = / 



Oi-y^)0{yo -y^) 



x{A\{p,{0) ) (yr ) V, (% ) I A) , (8) 

where ^ — yi — y2, yo, yi and y2 are space-time coordi- 
nates associated with the quark and gluon fields as illus- 
trated in Fig. [TJ The relative transverse coordinate is 
the Fourier conjugate of the transverse momentum qt in 
the gluon distribution function. 

Even though the above parton energy loss is derived 
for quark production and propagation in DIS, it is also 
valid for high-energy heavy-ion collisions. In the latter 
case, we assume the life-time of the quark-gluon plasma 
to be much longer than its formation time and expansion 
time scale and therefore can treat thermal partons inside 
the produced dense matter as in asymptotic states. One 
therefore can neglect correlation between the initial pro- 
duction rate of the jet parton and the quark and gluon 
density of the produced medium. The quark-gluon cor- 
relation function will then take a factorized form. 



dy' 



d^k 



(27r)32fc+ 



xe 



-.k+z 



,f(k,y)drd'^T 



dy' 



d^k 
(2^ 



fik,y) [l - cos{xLk+y )] 



X [4)k(xT + XL,qT) + c{xL)4'k{xT,qT)] , (9) 

where f{k,y) is the local phase space distribution of 
the color sources in the medium and c{xl) — fq{x + 
Xh)l jq{x) is the relative initial quark distributions in DIS 
and is given by the corresponding ratio of jet production 
cross sections in heavy-ion collisions. The unintegrated 
gluon distribution function per color source 4>kix,qT) is 
defined as 



4>kix,qT) 



27rfc+ 

x{k\F-+{0)F+iaT)\k). (10) 



The structure of the quark-gluon correlation function 
in Eq. ([9]) corresponds to two different bremsstrahlung 
processes and their interference associated with the 
different pole structures in Fig. [T] One can also catego- 
rize them according to how the longitudinal momentum 
transfer xl is provided. In the first term, the final gluon 
is induced by the secondary scattering with the medium 
gluon in which the intermediate gluon is off-shell. The 
longitudinal momentum transfer is therefore provided 



by the medium gluon and the contribution is propor- 
tional to gluon distribution 4>k{xT + x^^qr) per medium 
or "constituent" parton. These secondary processes cor- 
respond to quark-gluon Compton scattering where the 
initial gluon comes from a thermal constituent parton 
with a gluon distribution 4>k{xT + XL,qT)- Among these 
processes, one can identify a special case in which the 
quark scatters directly with a medium constituent gluon 
{xl = 1) as purely elastic processes and the correspond- 
ing energy loss as the conventional elastic energy loss 
[2^ . The second term in Eq. © corresponds to the pro- 
cesses in which gluon radiation is induced by the hard 
scattering that produces the initial jet (the quark after 
the photon interaction in Fig. [T]is off-shell) and the final 
gluon scatters again with a soft medium gluon with mo- 
mentum fraction xt- It is therefore proportional to the 
soft gluon distribution (j)k{xT, qr)- The longitudinal mo- 
mentum transfer x^ of the bremsstrahlung in this case is 
provided by the initial hard process with the cross section 
given in c{xl)- 

One can now define a generalized jet transport param- 
eter, 



qR(E,XL,y) 



A-K^asCR f d^k 



(2^)3 /(fc-y) 
(j)k{xT + XL^qr), (11) 



which includes the extra longitudinal momentum transfer 
XL from the medium to the propagating parton and the 
radiated gluon. The total parton energy loss from Eq. ([5]) 
can be expressed as 



dy dzdil^^^ ^ 



~c(xL)(7fl(£',0,y)]sin 



4y 



2^[qB.{E,XL,y) 

, (12) 



4Ez{l - z) 



in terms of the generalized jet transport parameter. The 
first term with the generalized transport parameter in- 
volves energy transfer between the propagating parton 
and the medium. It contains (but not limited to) what 
is normally defined as pure elastic energy loss [2a |. The 
second term that is proportional to the normal (or spe- 
cial) transport parameter qR^E^y) — qR{E,XL — 0, y) 
corresponds to pure radiative energy loss. 

Completing the integration over the phase-space of the 
radiated gluon, one can recover from the second term a 
similar form of total radiative energy loss in a static and 
uniform medium with finite length as in Eq. ([T]). How- 
ever, one needs to know the XL-dependence of the unin- 
tegrated gluon distribution function in order to calculate 
the "elastic" part of the energy loss. Furthermore, the 
transport parameter as defined in Eq. ([2]) should have 
some non-trivial jet energy (E) and temperature (T) de- 
pendence. 

Within the framework of twist expansion, the trans- 
verse momentum broadening of the quark jet has also 
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been calculated [2€ 



m - 1 



mx) 



= J dy-qn{E,0,y), (13) 



which is directly related to the normal transport pa- 
rameter qii{E,y) = qji{E,0,y) as defined in Eq. (fTTI) 
in terms of the unintegrated gluon distribution density 
of the medium. The jet transport parameter qii{E,y) 
therefore can be interpreted as the transverse momen- 
tum broadening per unit length for the propagating par- 
ton, as defined in Eq. Resummation of higher twist 
contributions leads to a diffusion equation for the trans- 
verse momentum distribution in which the above is the 
averaged transverse momentum broadening [27j . 

The approach leading to the above total energy loss 
and transverse momentum broadening has gone be- 
yond the conventional static potential model in two as- 
pects: (a) The result includes the longitudinal momen- 
tum transfer (xt) between the jet parton and the medium 
parton, which is related to the transverse momentum 
transfer (qx) through the unintegrated gluon distribu- 
tion density of the medium. This will result in the jet 
energy dependence of both the generalized and the nor- 
mal transport parameter qR,{E, y) (or transverse momen- 
tum broadening) which is absent in the static potential 
model, (b) The formula also includes the processes in 
which longitudinal momentum transfers xl comes from 
the medium parton and therefore it depends on the gen- 
eralized jet transport parameter qnixL)- It contains con- 
tributions from elastic energy loss. 

For the remainder of this paper we will focus on the 
energy dependence of the normal transport parameter 
qji = qii{E,y). Since it is essentially the transverse mo- 
mentum broadening per unit length which can be directly 
measured in experiments such as DIS and 7-jet events in 
heavy-ion collisions, wc will suppress the space and time 
dependence to simplify the notation. 



III. GLUON DISTRIBUTION IN A 
QUARK-GLUON PLASMA 

As shown in Eq. (fTTj) . the transport parameter qu 
experienced by a propagating parton can be defined in 
terms of the unintegrated gluon distributions 4>k{x,q^) 
of the color sources in the quark-gluon plasma. After av- 
eraging over the momentum of the color sources, it can 
be expressed as, 

qR = — tP I ~o / dx 



x6{x - 



2p' (fc^ 



-)a,(g|)0(x,(z^), (14) 



where is the average energy of the color sources and 
(f>{x, q^) is the corresponding average unintegrated gluon 



distribution function per color source. The integrated 
gluon distribution is 



xG{x, fj,^) 



d^qr 
(27r)2 



[x^qr). 



(15) 



We have extended our earlier definition of qa to include 
the case of a running strong coupling constant in 
QCD. We will refer to the case of fixed coupling constant 
as conformal gauge theory. However, for any scale below 
the temperature /i^ < we will consider frozen and 
treat it as a constant. 

Consider the lowest order parton-parton small angle 
scattering, we can obtain q^ as, 



d^k 

d^k' d^p' 



^f,[k){l± h{k'))ql\Mm? 
(27r)V(p + fc-p'-fc'), (16) 



where AAm is the truncated parton-parton scattering 
matrix element, 



M 



Rb 



1 



(1 — Xg) C0S( 



where cos ( 



f (1 - X^) + flj^lTTiXq) + ^^ag 

[v X q) ■ {vb X q)lq^ , Xq 



(17) 



qo/q 



and fj,jj = g (Nc + n//2)T^/3 is the Dcbye screen- 
ing mass in thermal QCD medium with temperature 
T and /imag ~ ^c5^/27r is the non-perturbative mag- 
netic screening mass [281 , [29l . [30| . The color factors 
for different scatterings are Cqq = Cf/2Nc, Cqg = 1/2, 
Cgg = N^/{N^ - 1). The statistical factor i^t, is 2{N^ - 1) 
for gluons and ANcUf for n/ flavors of quarks. We use an 
effective gluon propagator to include the resummation of 
hard thermal loops (HTL) [3^. The scaled self-energies 
in the effecti ve p ropagator in the long-wavelength limit 
are given by [3^ . 



TTLiXq) = l-:iiln'^^ 



1 +2^0 \ , .71" 

+ l-Xq, (18) 



2 \1-Xqj 2 
-l-Xq{l - Xg) . 



(19) 



One can rewrite the phase-space integration in Eq (jl6[) 
as 

d^k' d^p' 4 



(27r)3 (27r)3 



{2nyS\p + k - p' ~ k') 



(27r) 



^ dxd^qrSix- -^), (20) 



2p-k+- 



where x = q^/k~^. For small angle scattering, one can 
set q^ w q^ and Xq ~ xk~^/qT. We further approximate 
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fc+ by its average value = 3T in the scattering ma- 
trix. Note that energy-momentum conservation fixes the 
relative angle between k and q. Therefore, the angular 
phase-space for k is only 27r. One can complete the rest of 
the phase-space integration over the initial momentum, 

^/,(fc)(l±/,(fc')) - j^h{k){l±h{k)) 

= (gluons) or — (quarks). (21) 



Using 



t;Cr,i^, + CngVg = 2N,Cr{\ + ^) , (22) 



and 

P = 

one can express Eq. pB]) as 



K^m^9+^-^) - 2{Nl-l){l+:^) — C{-i), (23) 



-1 



9fl 



iV2- 1 



1 + nf/2Nc 



2tt 6C(3) l + 3n//4C; 



X / dxdq^S{x 



2p- 



-)ql\M 



Rb\ 



(24) 



The factor tt'^ /6({3) comes from the quantum statistics 
effect for the final state partons in the scattering pro- 
cesses. According to the definition in Eq. (fH|) . one can 
obtain the unintegrated gluon distribution function. 



1 + nf/2Nc 



6C(3) 1-K3n//4CF 



IMml'q^, (25) 



in a quark-gluon plasma and the integrated gluon distri- 
bution function is 



NcOs TT^ l + nf/2Nc 
27r 6C(3) 1-H3?i//4Cf 

X / dq^lMRbl^qr- 
Jo 



(26) 



We concentrate in the small Xq region, Xq — SxT/qx <C 
1. For a typical momentum transfer of the order of 
the requirement that Xq is small leads to a; ^ 
(\/^c5^/3)/3. Within this approximation, we obtain 

TTLiXq) ~ I - IX- , TTT(Xq) M -IX- , (27) 

Zqx ^qx 



and 



(x,g|,) = ..2 ^^,oN -, , o„ /.^ Hx,yq); (28) 



Ail^ 6C(3) H-3n//4Ci= 



1 



297r2r2/4^2^ 



(29) 



where = q'^/fijj. For a; 4/Lti„ag/7r = 2Ncg'^ /-k'^, one 
can neglect the magnetic mass and complete the integra- 
tion in Eq. ([26|) and obtain 

^2 



xG{x,^i') « iV, 



as TT^ l + n//2iVc 
2^6C(3) l + 37i//4Ci. 



ln(l + -^)^ 

3xr 

1 - 0.035 



D 



6 



ln(l 



Pd 
16 



1 + mVMd 

-3xTfj.o/l^''' 



P" 



9tt 



2 a;2T2^4^ ^ 



(30) 



where the first term is an approximation of the numer- 
ical integration from the electric part of the interac- 
tion for X < 2/X£)/37rT — 2y/ Ncg"^ /3/3tt. Because of 
the static Debye screening, it has a very weak x de- 
pendence in this X region, which can be ignored for 
large values of fi/fJ-o > 1- The magnetic part of inter- 
action, on the other hand, has only dynamical screen- 
ing and therefore lead to the dominant x dependence 
of the gluon distribution from a quark-gluon plasma at 
small X. However, our approximations are not valid 
for X < Afi^ag/'^T^T fa 2Nc Q^ / Sn"^, where the the non- 
perturbative magnetic mass ]28l. [29l [33| jimag ~ -^cff^/27r 
becomes important. In this region, the logarithmic a;- 
dependence of the gluon distribution from the magnetic 
interaction disappears and is replaced by a constant 

In(MVMmag)- 

For large ^^//^|) 3> 1 in the small Ncg'^ < x < y/Ncg 
region of a pure gluonic plasma {uf — 0), the gluon dis- 
tribution per gluonic color source is then. 



xG{x, /i^) « Ca- 



TT 6C(3) 2 



3i M +il,/f£ 
2 njj 3 xT 



(31) 



which is generated from perturbative gluon radiation. 
For a pure quark plasma, the corresponding gluon dis- 
tribution for each quark color source is 



xG{x, y?) f^ Gf' 



1 



TT 6C(3) 3 



o 2 

3 , M 

2 ^ll 



iln^ 

3 xT 



(32) 



For the remainder of this paper, we will focus on a pure 
gluonic plasma. 



IV. GLUON SATURATION IN A PLASMA 

Similar to gluon saturation in a large nucleus at small 
cc, saturation could also happen in the small a; region of 
a quark-gluon plasma. The saturation scale is given by 



Ql{x) 



47r2iVcQ!., 



-pxG{x, Q1) min(L, Lc) 



(33) 



where L^. — I /xT is the coherence length for parton scat- 
tering in a thermal medium. Since the HTL resumma- 
tion does not include coherence effects, the use of the 
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gluon distribution in Eq. (|3ip requires that the mean- 
free-path of thermal gluons must be larger than the co- 
herence length. Given the perturbative expression of the 
mean- free-path po| . 



this implies, 



977 " " Ncas 



A/ 



4C(3) N^al 
97r X 



In- 



Nra, 



or 



9tt 



< 1, 



Ncas 



(34) 



(35) 



(36) 



'-IX 

c 









oaiUialcO 












Linear 




:in(lsPg^T2) 


ln(N g^T^) 


lnT= 



In this regime, one can use the perturbative gluon distri- 
bution [Eq. (I3T|) ] to determine the saturation scale. 



pxG(x,Ql)Lc 



X 



Inf + ^ln^ 

9 xT 



Neglecting the logarithmic terms, one can get a simple 
expression for the saturation scale in the perturbative 
regime. 



\{x)lT^ 



(38) 



FIG. 2; Illustration of the hierarchy of the saturation scale 
Q^(x) below the hard scale fi^ = in a weak coupling glu- 
onic plasma. 



To obtain the transport parameter qn at scale /i^ < 

(37) 

due to interaction with the gluonic color sources via 
exchange of HTL gluons, one has to complete the integral 
in Eq. A simple integration in Eq. (|24p without 

considering effect of gluon saturation gives 



14 9 9 q 

10 /i^ 



(41) 



Since = - N^g'^T'^, we note the following 

hierarchy of the saturation scale in a perturbative gluonic 
plasma 



{N.g^fT^^^i-'Jovx^l 



Qlix) 



T ,for X - N^g - -j;^, 



T2 



(39) 



as illustrated in Fig. [21 

In the calculation of qu for interaction among thermal 
partons, the typical Xm = Ql/{s) — Q^/18T^. One can 
determine the saturation scale at Xm from Eq. ([57)1 ^, 



9 3 /I , 18T 
2 y TT 



(40) 



It is interesting to note that the gluon saturation scale 
for interaction among thermal partons coincides approxi- 
mately with the Debye screening mass. Therefore, resum- 
mation of HTL effectively provides some kind of mecha- 
nism for gluon saturation in a thermal gluon plasma. 



as obtained by a previous calculation of q^ with dynamic 
screening [sj] ■ One can also obtain the above result from 
the integrated gluon distribution 



- j:p^P [xGix,fi )\,^^2^/^2 ■ 



(42) 



In principle, one should take into account the effect of 
gluon saturation in evaluating the transport parameter 
in the region q^ < Q'^{x). In this regime, we can follow 
KNL model [35| and assume the saturated unintegrated 
gluon distribution as a constant in qx, 



(l){x,ql) = 



2Nra, 



{X, 



6C(3) 



ql<Ql; 



(43) 



ql > 



where (f>{x,yq) is given by Eq. l|29)) and the saturation 
scale Q^{x) is determined by Eq. ([57)1 . Using the above 
model for gluon distribution in the saturated regime in 
Eq. (|14p . one can evaluate the thermal parton transport 
parameter. The result. 



TTN^aiT^ In 



(44) 



^ Note that the typical Xm ~ NcQ^; thus the breakdown of the 
approximation in Eq. I|31|l due to magnetic mass coincides with 
the onset of non-linear effects 



with Ql^ = IST^, is nearly identical to Eq. (gTD. This is 
because the dominant contribution to the parton trans- 
port parameter comes from > Ql for large Q^J l^'^j^, > 
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1 and therefore the effect of gluon saturation is neghgi- 
ble in the calculation of transport parameter for thermal 
partons. 

The similarity between results in Eqs. and ([H]) 
is also an indication that saturation effect is already 
present in the unintegrated gluon distribution function 
(j){x,q'^) in Eqs. ([Ml) and due to HTL resumma- 
tion. One can clearly see this by analyzing the uninte- 
grated gluon distribution (j){x,q^) [Eqs. ((28)) and (|29)) ] at 
X = qT/{s)- For large ^ Mdj 4>{x,q^) ~ ^/qx- The 
electric contribution to (j){x,q'^){x — q^/{s)) reaches its 
peak value ^ NcOts/lJ^D at w ~ and vanishes 
at q^ = 0. Without the magnetic mass, the magnetic 
contribution to 4>{x,q^){x = q^/{s)), however, continues 
to increase at q^ < Q'i ^ Mij and reaches a finite value 
(j){x,q^){x = qT/{s)) ~ at q^ <^ fj.jj/T'^. However, 

contribution to qa from this region of limited phase space 
is sub-leading in the leading logarithmic approximation. 



evolution equation, 
d^xG{y,0 1 



(47) 



where, following Ref. (36|, we have defined variables y 
and ^ as 2 



V = In - . 



(48) 
(49) 



The asymptotic solution to Eq. (|17|) leads to a growth of 
the gluon distribution function of the order exp(V2Cy), 
while resummation of terms in Eg. (|46|) leads to 
eyiY>{asNcy) and exp(^) respectively |36l |. Therefore, the 
DLA approximation is valid as long as 



EVOLUTION OF THE THERMAL GLUON 
DISTRIBUTION 



i « 
aNcy « 



(50) 
(51) 



The gluon distribution function in Eq. (jSlD was ob- 
tained via parton interaction in a thermal medium with 
a HTL resummed gluon propagator and, thus, is only 
valid for scales fj?' < . At larger scales, radiation 
of hard modes, i.e., partons with momentum k > T, is 
possible. These processes lead to the evolution of the 
gluon distribution which in vacuum is governed by the 
BFKL/DGLAP equations in the linearized regime. In 
the medium this evolution may be modified due to the 
interaction of the radiated gluons with thermal partons. 
However, since the medium effects are of the order of 
/iD << T, we neglect those at hard scales and use the 
vacuum evolution to determine the gluon distribution. 
The previous computation in Eq. (ISip serves as an initial 
condition of this evolution at fx^ = T^. 

Since we are interested in the determination of q^ at 
large jet energies, we need to know the unintegrated 
parton distribution (j){x,q^) in Eq. p4|) at small x ^ 
(qx) /QET. For a large path length, the typical total 
momentum transfer, qL, which will set the scale of the 
process, is also large. This is the regime of the double 
logarithmic approximation (DLA), in which the BFKL 
and DGLAP equations coincide [Sa . [STj . In this approx- 
imation, all terms enhanced by two large logarithms of 
the type 



asik^)NJn^\n- 
/i^ x 



(45) 



are resummed. Thus, the DLA approximation is valid if 
terms of the above type are larger than those of type 



asik^)N,ln- 

X 



asik^)N,\n^ 



(46) 



Note that the definition of ^ allows to describe simulta- 
neously the evolution of a conformal and non-conformal 
theory. For these two cases we have 



2N, , ln(QVA") 



for fixed a. 



for running 



(52) 



where the reduced t'Hooft coupling is A = 2Ncag/'K and 
6= {llN^-2Nf)/l2T:. 

The general solution of Eq. (|T7)) can be found by per- 
forming a Laplace transformation and is given by [stI . [ssj ] 



xG{x,Q^) 



2-Ei 



(53) 



where a is any real number larger than the real part of 
any poles of D{n). The corresponding Laplace trans- 
formation D{n) of the gluon distribution in Eq. (|3ip at 

g2 = ^2 ^ T2 (e - 0) is: 



D(n) 



dye-"yxG{x,T^ 



27r 6C(3)3 



-41n^ 



(54) 



where is the strong coupling constant as evaluated 
at a scale that is proportional to T^, since Eq. ([3T|) is 
obtained through scattering between thermal partons. 



The resummation of the terms in Eq. (|45p leads to the fixed) 



^ The definition we use is slightly diff'erent from that of |36|| and 
is more suitable for the description of a conformal plasma {as 
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For large values, the integral in Eq. ((53)) can be 
performed by saddle point approximation, yielding 



NcO^ TT^ 1 



27r 6C(3)3V^(2^y)i/4 



(55) 



The above evolved gluon distribution function grows 
rapidly (faster than a power) with the rapidity y. Thus, 
at large y non-linear effects become important leading to 
parton saturation. Similarly, we can determine the sat- 
uration scale Q^{x) from Eq. ([33ll with the above gluon 
distribution function xG{x,Q), 

Ql = B{x, Ql)imn{L, L,) exp { ^26^} , (56) 

where S,s = £.{Ql) and 



Bix,Q^,) 



1 7r3 N,as{Ql) 

9C(3) N^.-l V(26z/)^/' 

T2 



21n- 



1/2 



(57) 



In solving the self-consistent equation Eq. ([SH]) we will 
neglect the weak dependence of B{x,Qs) on x and 
and treat it as a constant as compared to the dependence 
in the exponent. This is an approximation we will take 
throughout this paper. 

We now can use the evolved gluon distribution function 
in Eq. (|55|) to compute the jet transport parameter as 
defined in Eq. (fT4|l . In the linear evolution region {q^ > 
Qg), the unintegrated parton distribution is computed 
by taking the derivative of Eq. ([55]) with respect to the 
scale. Keeping the leading term in ^y {i.e. consider only 
the ^y dependence in the exponent) we find 

d 

4>DLA{x,qT) = ^■K—YxG{x,q^) 
oqrp 

_ n 2/ as{q^)Nc , 2 



^ It 

Using Eq. dSS]), we find at q^ = Q^, 

(pDLA[X,Qs) 



7r2 pmin(i,Lc) V2$yy 



xG{x,q^). (58) 



(59) 



At scales q^ < Qj, Eq. ([55)1 is no longer valid since 
saturation effects take place which tame the growth of 
the gluon distribution function. Inspired by the KLN 
model of saturation [35j . we use a simplified model for 
the unintegrated gluon distribution function 



{x,qT) 



TT^ pmin(L,Lc) V2?s y ' 
<pDLA{x,qT) , 



q^>Ql. 



(60) 



We can then express qji in Eq. (|T4|) as 



1R 



Gb. 47rV 



dx 



(2^)' 
y.as{qT)Nc(l)DLA{x,Ql) + 



)2 

i max 



X S(x 



ql 



Ql 



(pgr 

(27r)2 



■)asiqT)Nc(t>DLA{x,qT) ,(61) 



where Q^ax ~ 6i?T. Integrating out the (5-function, we 
have 



^2 2 
Ini 

X 

2\n^^{xQl,,) 



dx 



min(L, Lc) 
Cr Att^p 

Nr. m - 1 



x/ dxNcasixQ'^^J(j)DLA{x,xQ 

max ), (62) 

where ^ Ql/Ql 



-12 /r)2 

imax ■ 



VI. CONFORMAL PLASMA 

We first examine the behavior of the saturation scale 
and jet transport parameter in a medium with fixed cou- 
pling constant. For a medium length L that is always 
larger than the coherence length for any jet energy, we 
find 



In 



Ql 



In- 



1 



(63) 



for small x. This means that both constraints in Eqs. ([50]) 
and ([51]) are fulfilled at small coupling A. We can then 
use the DLA approximation to describe the evolution of 
the gluon distribution function and evaluate the satura- 
tion scale and transport parameter at small Xm ~ 1/E. 
Note that the eikonal approximation is valid for distances 
such that the total momentum transferred to the probe 
qL <C Qmax^ sincc Qmax the momentum transfer for a 
large angle (90°) scattering. From this requirement and 
Eq. §7^ we find 

Ql 



XL < Lc 



)2 

■' max 



Qi 



(64) 



which is compatible with the weak coupling approxima- 



tion and L > Lc if Qmax ^ 



We determine the saturation scale by solving the self- 
consistent equation Eq. ([56]) . Treating as a constant 
and using the definition of ^ for fixed coupling constant 
[Eq. dSH)] at Ql{x^) and x^ = Ql{x^)/Q 



liax^ one ob- 



tains 



Qs {Xm) 




B 



2 + A Tn^ 



In- 



■j^^ ^ max 



9 

2 ^ max 



2A 



■In' 



B 
TV 



(65) 
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In the large energy limit, the above solution simplifies to 



(66) 



To compute the quenching parameter we study numeri- 
cally the integral in Eq. ([S^ and find that, for the infinite 
conformal plasma, it can be well approximated by 

Cr Ql{x^) 



Nc Tim\{L,Lc{x,-n)) TT 



-Us {QliXm)) 



X A^,ln- 



e(Q?(a^m))-fA^ca.(Q?(x™))ln^ 



.(67) 



This is a very good approximation for values of A > 1. 
For small A it approximates the exact integral within a 
factor 2 as long as Aln(Q^^^/T2) > 1. Substituting 
Eq. (HH) in Eq. ^ we find 



1 /A 



In- 



In- 



2 V 2 /in^in; 



(68) 



As expected [39|, the transport parameter is determined 
by the saturation scale. 

To determine the dependence on the coupling, we sub- 
stitute the definition of B and set — 



OR 
T3 




V2+I- VA 



4+ A 
2 + A 



36 

1 

4 



(69) 



Let us point out two interesting features in the above 
result: 1) qr grows as a coupling dependent power of the 
energy and 2) it depends non analytically on the reduced 
t'Hooft coupling A. The non-analyticity is a consequence 
of the evolution process. 

The derivation of Eq. (|69p for a conformal plasma is 
strictly valid for small values of A since both the evolution 
equation Eq. (|47)) and the initial conditions Eq. ((3T]) arc 
based on perturbation theory. However, in our computa- 
tion we have not make any further assumption about the 
smallness of A. Given the recent interest in the compu- 
tation of transport properties in strongly coupled Af = A 
SYM [U m, [ll, [iO, EH, it is stiU instructive to study 
the strong coupling behavior of the jet transport param- 
eter. Plotted in Fig. O is the jet transport parameter as 
a function of the reduced t'Hooft coupling A, normalized 
by its large coupling limit, 

qR[X = 00) = — = — — . (70) 



II 
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FIG. 3: Normalized jet transport parameter as a function of 
reduced t'Hooft coupling A in a formal plasma for an energetic 
probe of log.oiQlaJT^) = 90. 



The normalized jet transport parameter increases mono- 
tonically with the coupling A and reaches its asymptotic 
value in the strong coupling limit. Note also that the 
above limit assumes that the energy of the probe is large 
such that ln((5^„^/r2) >> in A. 

Several comments on this strong coupling limit are in 
order: 

• In the strong coupling limit, the saturation scale 
approaches its maximum limit Ql{xm) — Qmax 
and, thus, the eikonal approximation is question- 
able. Both the saturation scale and the transport 
parameter become independence of the initial con- 
dition as contained in B. 

• Eq. (j70p has a power dependence on the energy of 
the probe, and the power becomes coupling inde- 
pendent in the strong coupling limit. 



Nc 



• Contribution to Eq. ((70)1 comes completely from 
the saturated part of the gluon distribution func- 
tion. We have performed a simplified treatment of 
this region by considering it constant. This is well 
motivated by numerical solutions of the Balitsky- 
Kovchegov equations at weak coupling [13, l43| . 
However, at strong coupling extra dependencies on 
the coupling (subleading at weak coupling) may be- 
come important. 

• We have not considered the evolution of the wave 
function of the probe. This is motivated by the 
weak coupling picture, in which such evolution is 
considered separately as the radiative processes of 
the probe and are described by radiative energy 
loss. At strong coupling this separation of the 
probe and medium evolution becomes ambiguous 
and may lead to an extra coupling dependence. 
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VII. NON-CONFORMAL PLASMA 

From the analysis of a conformal plasma with a fixed 
coupling we concluded that the saturation scale grows 
faster than any logarithmic jet energy dependence. Since 
the typical momentum scale is dictated by , effects of a 
running coupling constant become important in the QCD 
plasma for large energy probes. This issue is addressed 
by solving numerically Eq. ([5^ with ^ given by Eq. ([5^ . 
In this case we find that Eq. ((62)) is well approximated 
(within 20 %) ^ by 



In^ 



min(L,Le) In ^ 



A2~ 



In ^ In 



(inS/m^) 



Infln^ 



A2 



ln(l/x„) 
ln(Q2(2,^)/A2) 



, (71) 



where (5l = 1/2 if L > and (5^ = 1 otherwise. As in 
the conformal case, qB. is determined by the saturation 
scale, which is given by 

Ql = i?(x,Q2)min(i,L,) 



X e 



- In 1 



(72) 

where B{x,Ql) is given in Eq. ([57| and Lc = l/x„iT = 

As in the conformal case, the saturation scale and jet 
transport parameter have a fast growth with the jet en- 
ergy. In spite of the fact that the above results are de- 
rived with an approximation for asymptotically small x 
(which implies large saturation scales), we would like to 
make some numerical evaluations of the jet transport pa- 
rameter for jet energies accessible at RHIC and LHC and 
address the experimental consequences of this growth. 

We solve numerically the self-consistent equation 
Eq. ([7^ for the saturation scale Ql{xm)- In order to 
avoid the infrared singularity of as we regulate the cou- 
pling constant as 



a.(g') 



1 



1 



(73) 



To be consistent, we also replace 



In In 



A2 



'In' 



A2 



^ For determining this expression we assumed that the coherence 
length is always larger or smaller than the path length. In the 
numerical computations presented this is not assumed and the 
min() is replaced by a smooth function. 
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FIG. 4: (Color online) Coherence length times temperature 
as a function of the energy of the probe for different temper- 
atures. 



The coupling constant at a thermal scale is determined 
by solving the coupled equations. 



1 



1 



QliT) = 



18. (^T^ 



(74) 
(75) 



which are essentially Eq. ^ ai x = Ql{T)/l8T'^ [we 



have set the logarithms in Eq. ([37)1 to be of order 1]. 
Finally, since Eq. ([72]) is only valid for asymptotically 
large rapidities y, we also shift the rapidities to 



y ^y + VQ, 



(76) 



with j/o — 0.24. This value has been chosen such that 
as y decreases we recover the value of the saturation 
scale Ql{T). In the following numerical evaluation we 
wiU choose /i^ = and A = 200 MeV. 

Since the medium is finite in heavy ion collisions, we 
start by studying the coherence length. This is computed 
by solving 

1 



y2 



6E 



■In 1 



hi ■ 



(77) 



This coherent length is used to calculate the saturation 
scale for any large medium size L > Lc- For small 
medium size, L < Lc, the actual length L is used to 
calculate the saturation scale. The coherence length is 
plotted in Fig. |4] and it shows a strong energy depen- 
dence, as can be inferred from Eq. (j77p . This strong jet 
energy dependence Lc ^ VE is approximately indepen- 
dent of the evolution of the gluon distribution function, 
and stems from the definition of the coherence length 
as Lc = l/xT. As expected from the running coupling, 
it does not scale with temperature. For a characteristic 
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FIG. 5: (Color online) Saturation scale as a function of the 
jet energy. 



FIG. 7: (Color online) Saturation scale as a function of the 
path length for different probe energies. 
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FIG. 6: (Color online) Jet quenching parameter g as a func- 
tion of the jet energy. The square (triangle) marks the the 
value of q for thermal particle at T = 0.4 GeV [T = 0.6 GeV). 
Significant corrections to the energy dependence are expected 
at low energy which should approach their thermal value at 
E = 3T. 



temperature of T = 0.4 GcV in rclativistic heavy ion col- 
lisions, the coherence length is significant: for a probe 
of = 20 GeV, « 2.5 fm while at E = 100 GeV, 
Lc ~ 4.5 fm, which are comparable with the nuclear size. 

When the coherence length becomes comparable to the 
medium size a non-trivial length dependence of the sat- 
uration scale will arise, since the definition of is dif- 
ferent for path lengths longer or shorter than the co- 
herence length. This is illustrated in Fig. [5] where the 
saturation scale is plotted as a function of the energy 
of the probe. When the path length is longer than the 
coherence length, shows a stronger dependence on 
the energy. This is, in fact, driven by the energy de- 
pendence of the coherence length, and is mostly inde- 
pendent of the evolution of the medium gluon distribu- 
tion. When the path length is smaller than the coherence 
length we obtain a significant reduction of the saturation 



FIG. 8; (Color online) Jet quenching parameter q as & func- 
tion of the path length for different probe energies. 

scale and a much weaker dependence on the energy, since 
the DLA evolution leads to a growth that is weaker than 
a power but faster than a logarithmic dependence. We 
note that the gluon saturation scale obtained here for a 
gluonic plasma is significantly larger than that in a nu- 
cleon, where Ql ^ 1 GeV^ at a; « lO"'' 0. This is a 
consequence of the fact that the QGP is a much denser 
system than a nucleon or cold nucleus. The saturation 
scale in a heavy nucleus is enhanced by a factor of A^^^ 
and therefore might be large enough to facilitate a pertur- 
bative calculation of gluon distributions [i^. However, 
it is still an order of magnitude smaller than in a high 
temperature quark-gluon plasma. 

In Fig. [6] we show the value of the jet quenching pa- 
rameter qii from the integration of Eq. (|62p . For long 
path lengths q becomes path length independent. The 
leading energy dependencies of both Ql and Lc cancel 
and the observed energy dependence is a consequence 
of the evolution of the medium gluon distribution. At 
shorter path lengths, we obtain an enhancement of qn as 
a consequence of the evolution. 
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The length dependence of both the saturation scale 
and the transport parameter is shown in Figs. [7] and [5] 
for different energies of the probe. As shown in Fig.[7l the 
saturation scale grows as a function of the path length. 
However this growth is smaller than linear. Thus, the 
transport parameter, as shown in Fig. [51 diverges at 
small path length. Note, however, that at very short 
path lengths (L ^ A/) the mean momentum broadening 
should vanish, since the probe has no medium to scatter 
with. Thus, we expect correction to the small L depen- 
dence of both Ql and qR. Finally, when the path length 
is larger than the coherence length both quantities be- 
come length independent. Note that we have assumed a 
simplified transition from the region L < Lc to L > Lc- 
This is the reason for the abrupt change in the length 
dependence at L = Lc in Figs. [7] and [5] 

Let us remark that for path lengths smaller than the 
coherence length, the interaction of the probe with the 
whole length of the medium is coherent. Thus, if the 
length scales of space and time variation are smaller or 
comparable to the path length, the analysis of the satu- 
ration scales and jet transport parameter should be revis- 
ited. This will complicate the phenomenological extrac- 
tion of the transport parameter in an expanding medium 
with strong spatial variation as in semi-peripheral heavy- 
ion collisions. 

Because of the running of the coupling constant or the 
intrinsic scale (A) in QCD as a non-conformal gauge the- 
ory, the transport parameter qn has a non-trivial tem- 
perature dependence. To illustrate this, we plot in Fig. [9] 
the value of q^ scaled by the energy density, 



Stt^ 
15 



(78) 



to the power of 3/4 for the long path lengths (L > Lc). 
The dependence on the temperature is quite strong in the 
temperature range showed. This is not surprising since 
these temperatures are of the order of A and the coupling 
constant is very sensitive to the scale. The dependence, 
however, becomes weaker at higher temperatures. The 
jet energy dependence of the transport parameter is also 
stronger at lower temperatures. This, of course, is only a 
lower limit, since qa is larger for shorter path lengths 
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FIG. 9: (Color online) Jet quenching parameter scaled by 
e"^''*, with e the energy density. 



probe, 



A 



4p 



f 



(79) 



where we have used the definition of jet transport pa- 
rameter in Eq. ^ and (s) = Qmax ~ (^ET is the average 
center-of-mass energy squared of the jet-gluon scatter- 
ing. The requirement that the mean-free-path of the hard 
probe must be larger than the de Broglie wave length 1/E 
will set an upper bound for the energy loss parameter. 



QRiE) < 



's)E 3E^T 



2C 



(80) 



where C is a constant on the order of 0(1). 

We have checked that our numerical solutions of qn 
indeed satisfy this condition as shown in Fig. [TOl For a 
conformal plasma, q^ in Eq. (|69p increases monotonically 
with the reduced t'Hooft coupling A(see Fig. [3]). In the 
limit A ^ oo, Q'^{xm) = Qmax ^-iid the jet transport pa- 
rameter qji = CrTQ-^^^IANc for a gluon jet satisfies the 
bound for E > T. Since the strong coupling limit is an 
asymptotic behavior, the above bound on the transport 
parameter is also satisfied in the weak coupling limit of 
a conformal plasma. 

For thermal partons with (E) ~ 3T, the upper bound 
in Eq. ((80)) becomes 



VIII. BOUND ON q AND SHEAR VISCOSITY 
TO ENTROPY DENSITY RATIO 



Following Ref. [20|, one can relate the jet quenching 
parameter qn to the transport mean-free-path of the hard 



* The initial value of ijn before evolution is about half that of |46|| . 
The main reason is that we use as = 0.3 for T = 0.4 GeV while 
in Ei, as = 0.5. 



2 

> 



qniT) - 27C 



(81) 



According to Ref. [20'| , one can also relate the shear vis- 
cosity to the transport parameter for a thermal parton, 



77 « -sTXf 



3T3 



2qR{T)' 



or — ~ 

s 2qn{T) 



(82) 



Therefore, the upper bound on the transport parameter 
qniT) also provides a lower bound on the shear viscosity 



13 



m 1 

X 

a 10 

^ 10 



-2 



10 
10" 



-3 



10 



-5 




20 40 



60 



80 100 
E (GeV) 



FIG. 10: (Color online) Ratio iq/Ql^^^E = 2q/3E'^T as a 
function of jet energy E. 



to entropy density ratio 



s - 9C 



(83) 



This is very similar to the bound found by Danielewicz 
and Gyulassy from transport theory and the bound 
l/^TT found in the strong coupUng limit of A/" = 4 SYM 

The upper bound on the transport parameter qn in 
Eq. ([501) ^iid its connection with the shear viscosity in 
Eq. ([521) ^I'fi quite general since they do not rely on any 
assumption about the nature (perturbative or nonpertur- 
bative) of parton interaction inside the medium. They do 
rely, however, on a transport description of the plasma 
in terms of quasi-particles. Therefore, it is not surpris- 
ing that the connection between the transport parameter 
and shear viscosity does not hold in the strong coupling 
limit of A/" = 4 SYM theory, since thermal modes in such 
a strongly coupled theory cannot be described as quasi- 
particles [49| . 

We should emphasize that the jet transport parame- 
ter is an intrinsic property of the medium which could 
be dominated by non-perturbative physics. However, in 
the case of large saturation scale Qs ^"^d transverse mo- 
mentum transfer, the evolution of the gluon distribution 
function and the jet transport parameter should still be 
described by perturbative QCD and so should the in- 
teraction between jet and the medium and the radiative 
energy loss. Therefore, as far as the transport descrip- 
tion of the dense medium is valid, one can use the en- 
ergy dependence of the jet (parton) transport parameter 
Qr^E) as determined by jet quenching phenomenology to 
estimate the shear viscosity to entropy density ratio via 
extrapolation. 

With a recent phenomenological study of both sin- 
gle and dihadron spectra suppression using the next- 
to-leading order (NLO) pQCD parton model calculation 
[l^, the average gluon jet energy loss per unit length in 



a 1-d expanding medium is estimated to be 
dE 

1.9-3.4 GeV/fm, 



(84) 



for E = 10—15 GeV, which also includes an empiri- 
cal variation with jet energy. Using the relationship be- 
tween parton energy loss and the transport parameter in 
Eq. ([l]), one obtains an estimate of the average gluon jet 
transport parameter 



qoiE) 



ToasNc dL 



(^)irf « 1.0 - 1.9 GeVVfm (85) 



at an initial time tq — I fm/c Here, we used as ~ 0.24 
at ~ 5 GeV^ in a pure gluonic plasma. This roughly 
agrees with the numerical calculation in Fig. [51 As shown 
in Fig. [51 the energy dependence of qu is very weak 
for E < 20 GeV and long propagation length. Even 
for short propagation length L Lc ^ 2 fm, the en- 
ergy dependence is limited to about 25% variation. We 
therefore can use the above estimate for thermal parton 
transport parameter. Using the same initial temperature 
To = (337 ± 10) MeV as in Ref. one obtains 



m 

So 



3 T3 



2go(r) 



0.15-0.24. 



(86) 



For more consistent analysis, one should consider explicit 
energy dependence of the parton energy loss beyond that 
in qR. 

Finally, the requirement that the coherence length is 
larger than the mean-free-path of thermal particle sets a 
limit on the coupling constant a^. Indeed, from Eq. (|79p 



and Eq. dZH) qR 



2Qiixm)T, we find 



Xm = Qs{Xm)/{s) < 1/4. 

This bound is only satisfied at weak coupling 



(87) 



2Nr: 



27r In- 



Ncg^ 



< 1 



For larger values of the coupling , coherence effects 
in the multiple scattering of thermal particles become 
important. 



IX. SUMMARY AND DISCUSSIONS 

In this paper we have studied the energy dependence 
of the jet transport (or quenching) parameter qR. By re- 
lating (jR to the unintegrated gluon distribution function 



^ This phenomenological analysis was based on the assumption 
that q is independent of the path length. The nontrivial length 
dependence of q as obtained in this paper will affect the extracted 
average value of the jet transport parameter. 
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of the plasma we have shown that the energy dependence 
of qji arises from the evolution of the gluon distribution 
function. Thermal quarks and gluons in the plasma play 
the same role as the valence quarks of the nucleus and 
high energy jets probe their wave functions at small x. 
Similar as in a cold nuclear matter, the evolution leads to 
a growth in the gluon number which is eventually tamed 
by saturation effects. Therefore, the jet transport pa- 
rameter, also defined as the momentum broadening per 
unit length, is determined by the saturation scale Q^. 

Using thermal field theory with HTL resummation, we 
have derived the gluon distribution function for scales 
/Lt^ < as probed by the interactions among thermal 
partons. For such interaction among thermal partons, 
the coherence length is smaller than the mean-free-path 
and, thus, the saturation scale grows fast for small x, 
Ql ^ 1/x. Remarkably, evaluating the saturation scale 
at the typical x of the scattering among thermal partons, 
X « QI/AT'^, leads to ~ A*i)- Therefore, the typical 
momentum transfer is of the order oi as expected. 
What is more interesting is that, for large angle scattering 
of thermal particles with a; ^ 1, the saturation scales is 
of the order of the magnetic mass Ql ^ fJ-mag- 

The hard thermal loop result for gluon distribution 
function serves as the initial condition of the evolution 
of the gluon distribution as probed by an energetic jet. 
Since this process involves scales much larger than the 
medium scale, we have neglected thermal modification 
of the evolution kernel. We then used the double log- 
arithmic approximation to describe the evolution in a 
conformal theory and in (pure glue) QCD. 

For a conformal plasma, both the saturation scale 
and the jet transport parameter qn grow with energy 
as a coupling-dependent power. The evolution leads 
to a qn which is non-analytic in the t'Hooft coupling 
X — Nc- In the large coupling limit, qn becomes inde- 
pendent of the coupling and grows linearly with E. This 
is very different from results obtained in A/" = 4 SYM 
[U [11 HmO, SJ. As remarkable as this may be, the 
analysis we have performed is perturbative in nature and 
the extrapolation to infinite coupling might not be justi- 
fied. 

In the case of (pure glue) QCD, the evolution leads to a 
jet energy dependence of the transport parameter that is 
stronger than any power of logarithmic dependence. The 
saturation effect also gives rise to a non-trivial length de- 
pendence of the jet transport parameter. The running 
of coupling constant also results in a significant temper- 
ature dependence which becomes weaker at higher tem- 
peratures. We have numerically evaluated the saturation 
scale and jet transport parameter in a temperature range 
T = 0.4 — 0.6 GeV that is relevant for relativistic heavy 
ion collisions at RHIC and LHC. The growth of qa with 
jet energy is modest for large medium size L ^ Lc- How- 
ever, the energy dependence is significant for L ^ Lc- 
The obtained transport parameter is also larger than 
that computed via perturbation theory without evolu- 
tion [H]. For T = 0.4 GeV and = 10 - 20 GeV, our 



computed value qji « 1.5 — 2 GeV^/fm for a gluon jet 
is in agreement with the results from recent phenomeno- 
logical studies of experimental data on jet quenching at 
RHIC [iil, [H [H, Im^. It is, however, smaller than the 
results from phenomcnological studies [l^, [5l| based on 
an implementation of energy loss by Salgado and Wiede- 
mann (SW) j52j . A recent study within this model with 
explicit space-time dependent profiles of energy density 
from 2 and 3-D hydro calculations gives a transverse aver- 
aged qR ~ 4—5 GeV^/fm [11] at initial time tq = 1 fm/c. 
Inclusion of dihadron suppression in the phenomenolog- 
ical study has been shown [ll] to greatly improve the 
sensitivity of jet quenching to variation of qr. It is clear 
that inclusion of the energy and length dependence of 
the jet transport parameter will also influence the phe- 
nomcnological study of the jet quenching measurements. 

Given the relation between shear viscosity rj and trans- 
port parameter qR,{T) for a thermal parton as derived re- 
cently in Ref. _20j, the energy dependence of jet transport 
parameter determined from theoretical and phenomeno- 
logical studies can also be used to estimate the shear 
viscosity in the dense matter produced in heavy-ion col- 
lisions. This will unify high pT and low px aspects of 
heavy-ion collisions. The latter can characterize the col- 
lective behavior of the produced dense matter which is 
well described by relativistic hydrodynamics with a negli- 
gible shear viscosity 54] ■ We have also derived an upper 
bound on the transport parameter of high energy jets. 
This upper bound can lead to a lower bound on the shear 
viscosity to entropy density ratio which is consistent with 
other transport studies. 

The saturation scale in a glue plasma we obtained in 
this study is much larger than that in a cold nucleus or 
a nucleon at extremely small x. This results from the 
high gluon density in a plasma with coherence length 
comparable to the medium size. This is quite different 
from the analysis of Hadron-Electron Ring Accelerator 
(HERA) data which leads to a power dependence of Ql 
on a: [i^l 



= IGcV 



3 X 10" 



0.288 



(89) 



where the coherence length is much larger than the nu- 
cleon size. When the coherence length is comparable 
to the medium length, the saturation scale is not lin- 
ear with the path length. Thus, instead of using the 
phenomcnological expression Eq. (|89)) . we used the DLA 
approximation to estimate Q^. Numerical solutions to 
the Balitsky-Kovchegov (BK) equation with running cou- 
pling constant [42,, J^] show that the unintegrated gluon 
distribution is consistent with the DLA asymptotes and 
that the saturation scale behaves as ^ exp{Ay^}, 
similar to what is expected from the DLA asymptotes. 
Note, however, that these conclusions are for large values 
of the rapidity, y ~ 10, while in our case the typical x 
probed by high energy jets in the kinematic range rele- 
vant to heavy-ion collisions at RHIC and LHC is not very 
small, X ~ 0.1. The large values of x obtained imply that 
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significant corrections to the obtained behavior may oc- 
cur, which could be addressed by a numerical analysis of 
the BK equation. 
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